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Abstract 

Mathematical models related to some Josephson junctions are pointed out 
and attention is drawn to the solutions of certain initial boundary problems 
and to some of their estimates. In addition, results of rigorous analysis of the 
behaviour of these solutions when t —>■ oo and when the small parameter e 
tends to zero are cited. These analyses lead us to mention some of the open 
problems. 

Keywords: Third order parabolic operator; fundamental solution; supercon¬ 
ductivity; Josephson junction. 

82D55, 74K30, 35K35, 35E05. 

1 Introduction 

Our purpose is to: 

i) furnish a short review of the mathematical contributions to the dynamics of 
the Josephson junctions, 

ii) introduce some possible open problems. 

From the mathematical point of view, many descriptions of superconductivity 
phenomena have been developed and an important contribution has been given by 
Brian David Josephson. He predicted in 1962 the tunnelling of superconducting 
Cooper pairs through an insulating barrier to pass from one superconductor to 
another (Josephson effect). He also predicted the exact form of the current and 
voltage relations for the junction (Josephson junction) [1], (Experimental work 
proved that his theory was right, and Josephson was awarded the 1973 Nobel Prize 
in Physics.) 
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The flux-dynamics of a Josephson junction, i.e, two layers of superconductors 
separated by a very thin layer of insulating material, can be described by means of 
Sine Gordon equation (SGE): 

( 1 . 1 ) Uxx-utt = smu, 

where x denotes the direction of propagation, t is time and the variable u = u{x,t) 
represents the difference between the phases of the wave functions of the two super¬ 
conductors. 

However, in dealing with real junctions it seems necessary to take into account 
other effects such as losses and bias. Therefore, many authors prefer to consider the 
so-called perturbed sine Gordon equation (PSGE): 

(1.2) euxxt + uxx-utt-aut = sin u-'j. 

In this case, terms euxxt and aut represent respectively the dissipative normal 
electron current flow along and across the junction, (longitudinal and shunt losses) 
while 7 is the normalized current bias [2]. The value’s range for a and e depends 
on the real junction. Indeed, there are cases with 0 < a, e < 1 and, when the shunt 
resistance of the junction is low, the case a large with respect to 1 arises [2-4]. 

In some cases, extra terms must be considered. Eor example in a semiannular or 
in a S-shaped Josephson junction, when an applied magnetic field b parallel to the 
plane of the dielectric barrier is considered, the dynamic equation is; 

(1.3) £Uxxt + Uxx — utt — CLUt = sintt — 7 — bcos{kx), 

where the last term evaluates a transient force on the trapped fluxons and locates 
these ones at the center of the junction. [2,5,6]. Moreover, if an annular junction, also 
provided with a microshort, is considered, the vortex dynamics in a static magnetic 
field is modelled with the general perturbed sine-Gordon equation (see,f.i. [7]): 

(1.4) euxxt + Uxx - Utt - aut = [1 - S{x)fj] sin tt - 7 - bcos{kx), 

where /r is the current density associated with the microshort. 

Nowadays, in addition to rectangular or annular junctions, many other geome¬ 
tries for Josephson junctions have been proposed. For instance, window Josephson 
junctions (WJJ) ( [ 8 ] and reference therein) or exponentially shaped Josephson junc¬ 
tions (ESJJ). [9-12]. This type of junction is only a particular case of a structure 
covering a region 


0 < X < L, g 2 ix) < y < gi{x). 


0 < w{x) = gi{x) - g 2 {x) < 1 , 


(1.5) 

Denoting by 

( 1 . 6 ) 
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the evolution of the phase inside the junction is given by: 


(1.7) 


£Uxxt + Uxx - Utt - aut = sin u- r(x) - 


'w{x) ^ 
w{x) 


£Uxt) 


w{x) 

w[x) 


where T{x) = and rjx rjy is the normalized magnetic field respectively in 

the X and y directions [10]. When one assumes gi{x) = —g 2 {x) = Woe~^^, where 
A is a constant that, generally, is less than one, an ESJJ is obtained. Moreover, 
assuming that there is no bias current so that r(x) = 0 and gy = 0, the equation 
achieved is the following: 


( 1 . 8 ) 


£Uxxt + Uxx — Utt — £^Uxt — ^Ux — aut = sinu. 


The current due to the tapering is represented by terms A Ux and A e Uxt ■ In 
particular Xux characterizes the geometrical force driving the fluxons from the wide 
edge to the narrow edge. These junctions assure many advantages compared to 
rectangular ones, such as a voltage which is not chaotic anymore, but rather periodic 
excluding, in this way, some among the possible causes of large spectral width. It 
is also proved that the problem of trapped flux can be avoided (see f.i. [10]). 

There exist numerous applications of Josephson junctions especially as super¬ 
conducting quantum interference device (SQUID), which consists of a loop of super¬ 
conductor with one or more Josephson junctions. These devices are one of the most 
important applications of superconductivity. They are basically extremely sensitive 
sensors of magnetic flux. This peculiarity allows to diagnose heart and/or blood 
circuit problems using magnetocardiograms and even to evaluate magnetic fields 
generated by electric currents in the brain using magnetoencephalography -MEG- 
[2]. SQUIDs are also used in nondestructive testing as a convenient alternative to 
ultra sound or x-ray methods ( [2] and reference therein). In geophysics, instead, 
they are used as gradiometers [3] or as gravitational wave detectors ( [4] and refer¬ 
ence therein). SQUIDs play an important role in the study of the potential virtues 
of superconducting digital electronics, too [13]. 


2 Mathematical Models and Equivalences 

All equations previously considered have something in common. More precisely, if 
one denotes by C the following linear third order parabolic operator: 

(2.1) C = edxxt - dtt + dxx - adt, 

(1.1) -(1.4) and (1.8) can be expressed by means of the unique equations: 

(2.2) Cu = f{x,t,u). 

According to the meaning of /, numerous other examples of dissipative phenomena 
can be considered. For example, equation (2.2) arises in the motion of viscoelastic 
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fluids or solids {see [14-17] and references therein} and in the study of viscoelas¬ 
tic plates with memory, when the relaxation function is given by an exponential 
function [18]. It can also be employed in the analysis of phase-change problems for 
an extended heat conduction model [19,20]. In addition, equation (2.2) arises also 
in heat conduction at low temperature [15, 21] and in the propagation of localized 
magnetohydrodinamic models in plasma physics [ 22 ]. Still, it is possible to find 
others in [23-26]. 

Then, an equivalence between the third order equation (2.2), typical of Josephson 
junctions, and biological phenomena has been pointed out in [27]. Indeed, let us 
consider the FitzHugh-Nagumo system (FHN) [28,29]: 


(2.3) 


d 


u 




= s 


d: 


dt 

— = bu - I3v, 
ot 


— V — au+u‘^(a + l— u) 


(0 < o < 1 ), 


where u{x, t) represents a membrane potential of a nerve axon at distance x and 
time t, and v{x,t) is a recovery variable that models the transmembrane current. 

This reaction-diffusion model characterizes the theory of the propagation of nerve 
impulses, and the connection between a third order equation like ( 2 . 2 ) and the 
(FHN) system can be realized changing the first one into the second one under 
continuous parameter variations [27]. 

An equation that is able to model all these physical problems has been introduced 
in [30] and it is represented by the following parabolic integro-differential equation: 


(2.4) 


Cr u = ut — euxx + au + b e u{x, r) dr 

Jo 


F{x, t, u). 


Indeed, it has been proved that (2.4) characterizes both reaction diffusion models 
like the FitzHugh-Nagumo system and superconductive models [30-34]. 

In particular, perturbed sine-Gordon equation (1.2) can be obtained by (2.4) as 
soon as one assumes 


(2.5) 

and F is such that 


^ ^ R ^ 

a = a -, 0 =-, p = - 

e e e 


( 2 . 6 ) 


F{x,t,u) = — e e [ senu{x,T) — 7 ] dr. 

Jo 


Furthermore, the integro-differential equation (2.4) is able to describe the evo¬ 
lution inside an exponentially shaped Josephson junction, too. Indeed, as it has 
already been underlined in [ 12 ], assuming 

(2.7) ^ ^ b= a/3=^-6. 
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with 



( 2 . 8 ) 


/i= e-^"[sin(e"^/2 


from the integro differential equation (2.4) it follows: 

(2.9) £Uxxt Utt Uxx (q: 4“ ^ ^ 

Therefore, assuming ^ u = u, (2.9) turns into equation (1.8). 


Remark: In (2.4) the kernel e u{x, r) can be modified as physical situa¬ 


tions demand and in this way many other physical phenomena could be described 
{see, f.i. [35-38]and references therein }. The particular choice made here is due to 
describe the superconductive and biological models considered. 

3 Mathematical Results 

There exist many significant analytic results concerning the qualitative analysis of 
equations related to Josephson junctions and many initial-boundary problems have 
been discussed in a lot of papers (see [15,39-43] and references therein). 

A first analysis, where the fundamental solution is determined, concerns operator 
C in case a = 0 [14,44]. Later, in [45,46], the fundamental solution of the whole 
operator C, of (2.1) is explicitly determined and various properties are analyzed. 
Estimates and properties of continuous dependence for the solution of initial value 
problem are determined, too. 

Moreover, in [47], in order to deduce an exhaustive asymptotic analysis, the 
Green function of the linear operator £ of (2.1) has been determined by Fourier 
series and by means of its properties, an exponential decrease of solution related to 
the Dirichlet problem is deduced. And still by means of Fourier series, existence and 
uniqueness for Dirichlet, Neumann and pseudoperiodic initial-boundary conditions 
are achieved,too [42,43]. 

The Dirichlet problem is still considered with respect to equation (1.8) and in [11] 
the problem is reduced to an integral equation with kernel G endowed with rapid 
convergence and exponentially vanishing as t tends to infinity. Indeed, let 



and 


(3.2) 


Gnit) 


1 


e sinh{unt), 
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the Green function is given by 


2 A ^ 

(3.3) G{x,t,C) = - Gn{t) siri'jnC sin^nX. 

n=l 

The initial boundary problem with Dirichlet conditions is analyzed and an ap¬ 
propriate analysis implies results on the existence and uniqueness of the solution. 
That is, indicating by 

Qt ^ {{x,t) : 0 < X < L ] 0 < t < T}, 


the following initial boundary problem 


(3.4) 


{dxx - >^dx) (eut + u)-dt{ut + au) = F{x,t,u), {x,t)£Q.T, 

u{x,0) = ho{x), ut{x,0) = hi{x), xG[0 ,L], 
u{0,t) = gi{t), u{l,t) = g 2 {t), 0 < t < T. 


for <71 = 52 = 0 admits the following integral equation: 

y —M 

(3.5) u{x,t) = {dt + a+eXdx -edxx) ho{^)e ^ G{x,^,t)d^ 

Jo 


+ [ hi{^)e ^^G{x,^,t)d^+ [ dr f G{x,C,t-T)e 2^T(^,r,u(^,r))d^. 

Jo Jo Jo 

So, a priori estimates, continuous dependence and asymptotic behaviour of the 
solution, are deduced, too. 

When boundary data are non null, in order to achieve explicit estimates of 
boundary contributions related to the Dirichlet problem, equivalence between the 
equation describing the evolution inside an (ESJJ) and the integro-differential equa¬ 
tion (2.4) has been considered. Indeed, operator Cr of (2.4) has already been exten¬ 
sively examined in [30] and the fundamental solution K with many of its properties 
have been determined. 

More in detail, if a, 6, e, fJ are positive constants, r = \x\ j^ and Jniz) 
denotes the Bessel function of first kind and order n, let us consider the function 


—— 1 ri~ ft 

(3.6) K{r,t) = ^-==e~‘^^--J— ^-==e“^(‘“*')ji(2, y 6 ^(t^^)d 5 . 

2y/net 2 V tte Jq yt — y 

The following theorem has been proved: 


Theorem 3.1. The function K has the same basic properties of the fundamental 
solution of the heat equation, that is: K{x,t) G G°° for t > 0, x G 3ft. 

For fixed t > 0, K and its derivatives are exponentially vanishing as fast as 
|x| tends to infinity. 
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For any fixed 6 > 0, uniformly for all |a;| > 5, it results: 


(3.7) 


lim K(x,t) = 0. 
UO 


For t > 0, it is Cr K = t). 
Moreover, it results 


(3.8) 


\K{x,t)\ < 


-jeI 

e 4et 

2 y/iret 


[ e 


e-at _ g-/3t 

fi — a 


Previous estimates show, as well, that K exponentially decays to zero as t in¬ 
creases. These and other properties also allowed to prove in [12] numerous properties 
of the following function which is similar to theta functions: 

OO OO 

(3.9) 9{x,t) = K{x,t)+ '^^[K{x+ 2nL,t) + K{x— 2nL,t)] = K{x + 2nL,t). 

n=l n=—oo 

So that, as for problem (3.4), denoting by 

G{x,f,t) = 9{\x-f\,t) - 9{x + f,t) 


and 


_ A 

F{x, t,u) = e 2 


f e [sin u) — 7 ] dr — hi{x) e 

Jo 


it has been proved that the problem admits the following integro equation: 


(3.10) u{x,t) = [ G{x,f,t)e 
Jo 

F XL 

-2e / 9x{x,t-T)gi{T)dT + 2£ - L,t - T)e~~g 2 {T)dT. 

Jo Jo 

Besides, a priori estimates and asymptotic properties have proved that when t 
tends to inhnity, the effect due to the initial disturbances {ho, hi) is vanishing, 
while the effect of the non linear source is bounded for all t. Furthermore, for large 
t, the effects due to boundary disturbances 51,(72 are null or at least everywhere 
bounded. 

Indeed, if /iq = hi = 0 and F = 0, the following theorem holds: 


"ho{f)df+ f dr f G{x,Ft)F{FT,u{x,T))d^ 
Jo Jo 


Theorem 3.2. When t tends to infinity and data gi {i = 1,2) are two continuous 
functions convergent for large t, one has: 


(3.11) 


u 


9l,oo 


sinh (To {L — x) 
sinh (To L 


+ 92,00 


sinh (To X 
sinh (To L ’ 


where ao = ^ and gi^oo = hmt^oo 5 i, (* = 1,2). Otherwise, when gi € Li[0, 00 ] (i = 
1 ,2) too, the effeets determined by boundary disturbance vanish. 
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Another aspect frequently highlighted in many papers is that the linear third 
order operator £ is an example of wave operator perturbed by higher order viscous 
terms. The behaviour of solution of (2.2) when a = 0, has been analyzed in various 
applications of artificial viscosity method [48,49]. Moreover, in [50], when e is van¬ 
ishing, the interaction between diffusion effects and pure waves has been evaluated 
by means of slow time et and fast times t/e. These aspects are also analyzed in [16] 
referring to the strip problem for equation (2.2) with a linear source term /, while 
in the non linear case, the Neumann boundary problem has been discussed in [51]. 

Also equation (1.8) can be considered as a semilinear hyperbolic equation per¬ 
turbed by viscous terms described by higher-order derivatives with small diffusion 
coefficients e. In [52], the influence of the dissipative terms has been estimated 
proving that they are both bounded when e tends to zero and when time tends to 
infinity, giving a mathematical proof of what has been observed in [9]. 

As for explicit solutions, an extensive literature exists, and more recently, various 
classes of solutions for (SGE) have been determined.(see,f.i, [53,54]). Furthermore, 
when e = 0, some travelling-wave solutions for (1.2) have been obtained both for jqj 
not larger than 1 and for jyj > 1 [55,56]. Still when e = 0, some classes of explicit 
solutions have been determined for equation (1.8),too [52]. 

4 Open Problems 

In light of what has been stated until now, many open problems can be highlighted. 

It would be interesting, for example, to study equation (1.2) when interface 
conditions for the phase (and its normal gradient) are added, connecting, in this 
way, with the problems of window Josephson junctions (WJJ) when the influence 
of an external magnetic field must be considered [57]. Indeed, letting e = 0, (1.2) 
exactly recalls one of the equations usually considered for (WJJ). 

When, on the other hand, e is not vanishing, a viscous term, represented by the 
third order term, appears. So that, it would be interesting to give an estimate of 
the diffusive effects due to the e-term,too. 

Moreover, according to the analogy between superconductor equations and reaction- 
diffusion models, the Robin boundary problem would be considered in order to 
achieve results for many biological phenomena, too [58,59]. 

Besides, as for analysis on asymptotic effects due to the boundary perturbations 
related to equation (1.8), as it has been pointed out, the Dirichlet boundary problem 
has already been considered in [12]. So, the evaluation could be extended to other 
boundary problems, such as, for instance, Neumann and mixed ones. 

Of course, in order to achieve estimates for other more significant physical prob¬ 
lems, this analysis and many other estimates could be carried out for solution of 
equation (1.3) and for equations like (1.4) where the presence of a gap in the vacuum 
chamber is considered, too [41] 

The analysis conducted so far required that in (2.4) constants a, b, e, (5 were all 


positive. This can be valid if we look for an analogy with an (ESJJ), but excludes 
application of (2.4) to some other junctions. Therefore it would be interesting to 
extend the analysis of operator Cr for any value of a, b, e, f3. 

Finally a qualitative analysis of operators should be made in case e, a, A were 
not constant. 

5 Conclusion 

The state of the art proves that many signihcant analytic results concerning the 
qualitative analysis of equations related to Josephson junctions have been obtained 
and many initial-boundary problems have been discussed. However other many 
important open problems may be considered and solved. 
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